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Abstract 

We calculate the dephasing rate of an electron state in a pinched quantum 

dot, due to Coulomb interactions between the electron in the dot and electrons 

in a nearby voltage biased ballistic nanostructure. The dephasing is caused by 

nonequilibrium time fluctuations of the electron density in the nanostructure, 

which create random electric fields in the dot. As a result, the electron level in 

the dot fluctuates in time, and the coherent part of the resonant transmission 

through the dot is suppressed. 
PACS numbers: 73.23-b, 72.10-d 
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I. Introduction 



The dephasing of electron states in Quantum Dots (QD) was considered mainly in con- 
nection with weak localization phenomena, see experiments and theory [Q,|J. A differ- 
ent type of phenomenon in which dephasing is important is interference phenomenon in an 
Aharonov-Bohm ring ||. If a pinched QD is embedded in one of the arms of such a ring 
the transmission through this arm is supported by a resonant electron state in the QD. The 
dephasing of this state || suppresses the interference in the ring, and this can be observed 
as a decrease of of the oscillating part of the ring conductance M . 

The dephasing is due to electron-phonon or electron-electron interactions of the QD 
electrons with some "environment" , which can either be in equilibrium or driven out of it by 
external forces. In the experiment |7j the dephasing was due to the capacitive interaction 
of the QD with a voltage biased point contact (PC), and the amount of dephasing was 
dependent on the bias. In a situation like this one can separate the equilibrium dephasing, 
which depends only on the temperature of the environment, from an additional dephasing 
which is due to voltages applied to the environment. The theory concerning this experiment 
was given in ||||. 

In recent experiments || the nanostructure (NS), that was capacitively coupled to the 
QD, was a multiterminal 2DEG device in a quantizing magnetic field . We present in this 
paper a generalization of the theory given in || that takes into account the specific effects 
appearing due to the complicated geometry, and the chirality of the states in the NS (see 
Fig.l). A similar problem was addressed in [TIJ using a different approach, based on lumped 



mesoscopic circuit elements. Broadening of electron transitions in self-assembled QD's due 
to Coulomb interactions with the electrons in the wetting layer was considered in [|TT] . 
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II. Model 



We consider a QD with a single level eo, that is described by the Hamiltonian Hqd = 
€qC + c, where c + is an operator creating an electron in the QD state. The NS is a multiter- 
minal junction in the 2DEG described by the Hamiltonian 

H NS = J dr^+(r)#(r)^(r), H(r) = (-iV - ^(r))" + U(r), (1) 

where U(r) is the potential confining the 2DEG, A(r) is the vector potential of the external 
magnetic field (e < 0, h — 1), and \I/(r) is the electron field operator. The capacitive 
interaction between the QD and the NS, assumed to be weak, is 

H mt = c + c J drW(r)p(r), (2) 

where p(r) = \P + (r)\I/(r) is the electron density operator, and W(r) is a Coulomb interaction 
kernel. 

The perturbation Eq. (^) has a "dual" meaning. If one combines H^s + H int one can see 
that W(r) is the change of the confining potential U(r) due to one electron occupying the 
QD state (when < c + c >= 1), while combining Hqd + Hi nt one can see that / drW(r)p(r) is 
the change of the energy eo due to the Coulomb interaction of the electron in the QD with 
the electron density in the NS. 

III. Dephasing rate notion 

At low temperatures the dephasing is due to electron-electron interactions [JP2"] . To 
calculate the dephasing rate we use the method developed in 0. The electron density in 
the NS fluctuates in time and creates a fluctuating potential in the QD which brings about 
fluctuations of the energy level eo- These fluctuations are given by 

5e (t) = J drW(r)5p(r,t), (3) 

where 5p(r,t) = p(r,t) — (p(r)). As a result, the correlator (5e (t)5e (0)) is defined by the 
density-density correlator {5p(r, t)Sp(r f , 0)), while {5eo) = 0. 
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Consider resonant transmission through the QD for electron energies e close to eo- When 
the QD level does not fluctuate the transmission amplitude t(e) contains the Breit-Wigner 
factor — i/[(e— eo)+«T], where T is the width of the level due to the QD's connection with the 
leads. When the level fluctuates the transmission and reflection can be elastic and inelastic. 
In interference experiments only the elastic transmission is of importance and to obtain the 
elastic transmission amplitude (t(e)) one has to replace the Breit-Wigner factor by || 

/ dtexp[-Tt-$>(t) + i(e-e )t}, (4) 

J 

where 

W) = \ f 'dH f dt"(6e (t')6e (t")). (5) 
2 Jo Jo 

One can see that |(i(e))| < |t(e)| which means dephasing of the QD state, that is responsible 
for resonant transmission. The same can be understood from the dynamics of the QD state 
amplitude §, (c(t)+c(0)) = (c(0) + c(0)) exp[ie t -Ft- $(t)]. 

The level fluctuations are characterized by their amplitude ((<5eo) 2 ) 1//2 and by the corre- 
lation time t c . The amplitude is proportional to the strength of the capacitive coupling W, 
while the correlation time is independent of W and is determined by the correlation time of 
the density-density correlator. Hence, for weak enough coupling, one has ((£eo) 2 ) 1 / 2 T c <C 1, 
which corresponds to dynamical narrowing. In this case $(£) = jt and the integral Eq.(£|) 
reduces to — i/[(e — e ) + i(T + 7)]. Here 

7 = vrir(0), (6) 

with the following level oscillations' spectrum 

1 f+co 

K{uj) = — / dte wt (8e (t)Se (0)). (7) 

This result means that in the case of dynamical narrowing one can describe the dephasing 
by a dephasing time t v = 7 _1 . The dephasing rate can be estimated as 7 ~ ((Se ) 2 )r c , and 
is smaller than the amplitude of the level fluctuations ((Seo) 2 ) 1 ^ 2 - In the general case when 
((<5eo) 2 ) 1 / 2 r c > 1, the transmission probability |(t(e))| 2 is not a Lorenzian, and a dephasing 
time can not be defined. 
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IV. Dephasing rate calculation 



To calculate the correlator K(u) we represent the field operator in terms of scattering 
states (SS's) |n[ (see Appendix) 

*( r ) = / ^E a -( e )^n(e,r), (8) 

where a^ n (e) is an operator creating an incoming electron in channel n of terminal a, with 
energy e. Performing calculations similar to those in || we find 

K(u) = ^K^(u), (9) 

act 1 

where the contribution from terminals a and a 1 is 

K™\u) = l -J^-J^-Y.f^m- U(e')]\W an , a , n ,(e, e')\ 2 [S(e' -e-u) + 5(e' -e + u)). 

nn' 

(10) 

Here we used the fact that when the SS's are normalized to a unit of incoming flux one 
has (a^ n (e)a a / n /(e')) = 2Tr5(e — e')5an,a.'ri f a (e), where / a (e) is the Fermi distribution for 
energy e in terminal a. It is convenient to write it as /(e — 5fi a ), where 5fi a = fi a — fi and 
/(e) = (e( e ~^/ T + l^j is the Fermi distribution with some reference chemical potential /i. 
The matrix element entering Eq.([T0|) contains SS's, 

WWn'M') = J drW(r)x a n{e,r)*Xa'n'(e',r)- (H) 

The integration here is over the interaction area, i.e. over that part of the NS which is close 
enough to the QD and where W(r) is not small (see Fig.l). 

In what follows we consider the case when the voltages V a applied to all terminals are 
small. We choose /x to be the equilibrium chemical potential (when all V a = 0) and 5fi a = 
eV a . In this case the relevant energies in Eq.(|H]) correspond to the small energy window 
|e — /j | ^ max[T, eV], where electron exchange between terminals happens. We assume that 
within this energy window one can neglect the energy dependence of the scattering states 
and hence of the matrix elements Eq.(|ll|). As a result we have 
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K aa (u) 



\\W a , al \ 2 J^J ^f(e - eV a )[l - f(e' - eV a ,)\W - e - U ) + 5(e> - e + «)], 

(12) 

with an effective matrix element 

\W a ^\ 2 = Y,\ W «n,a'n>\ 2 - ( 13 ) 



nn' 



Shifting the integration variables in Eq . (|12|) by eV a and eV a i one can see that the diag- 
onal contributions K aa do not depend on the applied voltages V m and are equal to their 
equilibrium values at V a = 0, i.e. 

K aa (uj) = J-^\W a , a \ 2 F T (uj), (14) 



where F T (u) = uj coth(cj/2T) = 2u \N t (oj) + §1, with N t (lu) = (e^ T + l) . Note that 
F T {0) = IT and F {u) = 

For the nondiagonal contributions a ^ a' we find after a shift of the integration variables 

K^ a '\u) = K aa '(oo) + K a ' a (ui) = J-\W a , al \ 2 [F T (u + eV aa ,) + F t (uj - eV aa ,)\, (15) 

where V aa > = V a — V a > . 

Using Eq. (Q) one can find the dephasing rate as a sum over single terminals and pairs of 
terminals, 

7 = E7 (a) +E7 M (16) 

a a<a' 

It follows from Eq.(|l4]) that a single terminal contributes to the dephasing only if SS's 
emitted from this terminal reach the interaction area, and that this is always an equilibrium 
contribution, 

7 {a) = 7rK aa (0) = ^\W a , a \ 2 T. (17) 

A pair of terminals contribute to the dephasing only if SS's emitted from both terminals 
overlap in the interaction area, 



7 (W) = K (-') (0 ) = ±-\W a , a ,\ 2 F T (eV aa ,). (18) 

When both terminals are at the same voltage, the contribution of this pair to dephasing is 
an equilibrium one. 

Note that Ft(uj) contains the "zero point fluctuations", but they do not contribute to 
the equilibrium dephasing rate, given by K(oS) at u — 0, hence for zero temperature there 
is no equilibrium dephasing. 

If one is interested in nonequilibrium dephasing one has to look only at pairs of termi- 
nals which are at different voltages, and which send scattering states that overlap in the 
interaction region. The nonequilibrium contribution of such a pair is 

C-O _ M _ 7 (W)| y=o = ±\ Wa , a ,\ 2 [F T (eV aa/ ) - F T (0)). (19) 

For zero temperature it reduces to 

7 { v a ' ] \T=o = ^\W a , a ,\ 2 \eV aa/ \. (20) 

Using the dual property of the interaction between the QD and the NS, we consider now 
W(r) as a small variation of the confining potential U(r) due to an electron occupying the 
QD. As a result the scattering matrix of the NS is changed according to Eq.( [A15| ) from S 
to S + 5S. Using in addition Eq. ([A6| ) we can express the matrix elements Eq.(|ll|) in terms 
of the scattering matrix variation 

W aneta > n > e = % Si3 m , an (e)*8Sf3 mta > n >(e). (21) 

/3 m 

Pinching the QD to the Coulomb blockade regime one can change the number of electrons 
in the QD one by one and measure the variation of the conductance matrix due to an 
additional electron in the QD |13|]. Since Gf3 a = J2 mn \S/3 m , a n\ 2 — <5/3,a, this is a way to measure 
(in case of simple enough NS geometry) the variation 5S and the matrix elements W. This 
procedure was performed experimentally for the simplest NS, being a one channel point 
contact. 
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V. Dephasing versus current noise 



Dephasing is closely related to current noise since current fluctuations are related to 



charge density fluctuations by the continuity equation. The results obtained in |L4| for the 
current noise can be presented in the following form 



e 2 



(SMI*) = — 2 £ A#F T (eV P p). (22) 

Here the left hand side is the u = Fourier component of the current cross-correlator in 
terminals a and a', 

Aaa' _ ST^ \a ao 1 /no\ 

A3' ~~ ^Pmfl'm'^P'm'^rm l zc V 

mm' 

with 

Apm^'m' = Spa&P'a&mm' ~ ^an,l3m^cm,l3'm' , (24) 

n 

where the scattering matrix is at e = //. One can see that a single terminal (3 contributes 
to (5I a 5I a r) only if SS's emitted from this terminal reach both terminals a and a', and that 
this contribution, given by the term with j3' — (3 is always equilibrium. Pairs of terminals (5 
and P' contribute to (5I a 5I a ') only if SS's emitted from each of these terminals reach both 
terminals a and a'. This contribution given by terms with (3' ^ (3 contains a nonequilibrium 
part. These conditions are very similar to those in case of dephasing. 



VI. Examples and discussion 

We consider first a simple one channel 2-terminal device with a gate and a QD (see 
Fig. 2). The scattering matrix of the gate (in the absence of an electron in the QD) is 



S 







r t 




t r 





cos 9e ia sin 9e llS 



sin 9e l @ cos 6e ta 



a — (3 = re — (a — (3), 



(25) 



where r and t correspond to reflection and transmission of the SS's approaching the gate 
from left, while f and t correspond to the SS's approaching the gate from the right. If the 
magnetic field B ^ the scattering matrix is not symmetric, t ^ t. 



S 



In case of zero temperature the equilibrium part of the QD state dephasing rate vanishes, 
and the nonequilibrium part is according to Eq. (p0|) and Eq.(^TJ) 

1 = 7r \W\ 2 eV J \W\ 2 = \r*St + t*Sf\ 2 = \56 - i(5p - Sa) sin 6 cos 6\ 2 . (26) 

Here V = | V12I an d St , Sr are the changes of the transmission and the reflection amplitudes 
due to the electron in the QD. This result was obtained in |J for B = and a symmetric 
gate. The shot noise in this device (in the absence of an electron in the QD) is ]15| 



((Siy) = (e 2 /4n 2 )\r\ 2 \t\ 2 eV. (27) 

Both the dephasing and the shot-noise are due to the same nonequilibrium fluctuations, 
but they are not proportional to each other. To get some insight, consider first zero or a 
weak magnetic field, when both SS's 1 and 2 occupy the whole crossection of the sample 
(see Fig. 2). If we assume there is no reflection from the gate, i.e. \t\ 2 = 1, \r\ 2 = 0, we find 
((SI) 2 ) = 0, while 7 = 7r|5r| 2 eV 7^ 0. The dephasing is nonzero because SS's emitted from 
different terminals overlap near the QD, while the shot noise is zero because each terminal 
(where the shot noise is measured) is feeded only by one SS. The situation for the shot noise 
changes if the gate is reflecting, in which case each terminal is feeded by both SS's. 

One can understand this difference from the following simple calculation. In a channel 
without reflection the wave function is ip(x) = ae lkx + be~ tkx , where the two terms are SS's 
coming from the left and right terminals. The corresponding charge and current densities 
are p(x) = e{\a\ 2 + \b\ 2 + (ab*e l2kx + c.c.)} and j(x) = ev k {\a\ 2 — |6| 2 }. What is important 
for nonequilibrium fluctuations is the overlap of SS's coming from different terminals, i.e. 
terms proportional to ah. Such terms do not exist in j but do exist in p. This is why the 
shot-noise is zero, while the dephasing rate is not. The term (ab*e l2kx + c.c.) is of quantum 
origin. It means that in a quasiclassical situation, when the number of channels is large, 
this term will average out due to "integration" over k. When the gate is reflecting, \t\ 2 7^ 1, 
\r\ 2 7^ 0, both p(x) and j(x) contain terms proportional to ah. One can easily check it using, 
for example, the wave function to the left of the barrier ip(x) = a[e lkx + re~ tkx ] + bte~ lkx . 



It is also important to notice that for u = the current and charge fluctuations are not 
coupled by the continuity equation. 

Consider now the same device in a strong magnetic field when the SS's are edge states 
(ES's) localized near the boundaries. We assume also that the QD is far from the gate 
and the interaction region does not reach ES2. In this situation due to the chirality of ESI 
the QD can change only the phases of t and f. As a result \W\ 2 = (5j3) 2 \t\ 2 \r\ 2 , i.e. the 
dephasing rate is proportional to the shot-noise. This is because for chiral states the current 
and charge densities are proportional, j = pvk- (The connection between dephasing and the 
phase of t was mensioned in ||16|| ). 

As a second example we consider a 4-terminal device similar to that used in experiment 
||, with a geometry as shown in Fig. 3. The source S (a = 1) and the drain D (a = 2) are 
used to bias the device. Two floating terminals, one down-stream from S to D (a = 3) and 



one up-stream from D to S (a = 4) are "dephasors" according to [[[7]] . Gate A regulates the 
source-drain current, while gates B and C block the floating terminals. The QD is located 
far from gate B. We assume there is only one LL at the Fermi energy and that the ES's 
at opposite edges are well separated and do not overlap. We will be interested only in 
nonequilibrium dephasing and consider zero temperature. 

The SS's emitted from the up-stream floating terminal 4 do not reach the interaction 
region and hence this terminal does not contribute to the dephasing of the QD. (In what 
follows it is assumed that this terminal is blocked). Only scattering states emitted from 
terminals 1,2, and 3 overlap in the interaction region, and hence in accordance with Eq.(|T6D 
one has 7 = 7^ 12 - ) +7^ 23 ) + 7^ 31 ). Since the QD is located far from point contact B, all SS's in 
the interaction region have the form of the same ES if^(r) = w(r) with different amplitudes, 
i.e. Xi — e^HArsW, X2 = e^f^re-u;, % 3 = e l< ^H B w. Here r^, t^ and f^, Ia are the reflection 
and transmission amplitudes for ES's approaching A from left and from right, r^, t# and fs, 
ts correspond to ES's approaching B from above and below. The phase factors e*^ depend 
on the position of the QD. The relevant matrix elements Eq.(|I^) are 
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\W l2 \ 2 = \t A r B \ 2 \r A r B \ 2 \W\ 2 , \W 13 \ 2 = \t B \ 2 \t A r B \ 2 \W\ 2 , \W 23 \ 2 = \t B \ 2 \r A r B \ 2 \W\ 2 , (28) 
where 

W = J drW(r)\w(r)\ 2 . (29) 
Using these matrix elements we find 

7 {12) = A\t A \ 2 \r A \ 2 \r B \*\V 12 l 7 (23) = A\t B \ 2 \r B \ 2 \r A \ 2 \V 23 l 7 (31) = A\t B \ 2 \r B \ 2 \t A \ 2 \ V 13 \, (30) 

with the constant A = (e/47i)\W\ 2 . 

When terminal 3 is open, i.e. r B = 0, SSI and SS2 are absorbed in this terminal and 
then the interaction region is reached only by SS3. There is no overlap in the interaction 
region of SS's emitted from different terminals and as a result all the contributions to 
nonequilibrium dephasing rate vanish. When terminal 3 is blocked, i.e. t B = we find 

7 = 7 (12) =A |^|2| rA |2| ri2 |^ 7o . 

Since terminal 3 is floating V 3 is given by the condition that the current entering this 
terminal is zero, which leads to V 3 = Vi\t A \ 2 + V 2 \r A \ 2 . Using this one finds 

7 = 70 M 2 (2- |r B | 2 ). (31) 

One can see from this result that 7 < 70, i.e. the floating terminal suppresses the nonequi- 
librium dephasing rate of the QD state. This result is in agreement with experiment |J. We 
would like to stress that the suppression is not because of dephasing the SS's coming to the 
interaction region. The absolute values of the matrix elements that enter the expression of 
the dephasing rate according to Eq. (|20|) do not depend on the phases of the SS's overlapping 
in the interaction region. If one would simply destroy their phases it would not affect the 
dephasing rate 7. The floating terminal suppresses 7 because it absorbs the SS's moving 
towards the interaction region from different terminals. It is important to have in mind that 
a theory based on the representation Eq.(|8|) assumes that terminals absorb incoming waves 
as black bodies, which means that terminals have infinite capacitance. 
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It is instructive to compare the dephasing rate with the shot noise. When terminal 3 is 
blocked the shot noise is known to be ((Six) 2 ) = ((5I 2 ) 2 ) = (e 2 /47r 2 )e| Vu\ |*a| 2 |?"a| 2 - Using 
Eq. one can see that opening terminal 3 does not change ((5Ii) 2 ) but suppress ((5I2) 2 ) 
by exactly the same factor |t\b| 2 (2 — |re| 2 ) as the dephasing rate. 
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APPENDIX: 

In this Appendix we list some useful properties of the Green function, the scattering 
states and the scattering matrix, valid also when the magnetic field B ^ 0. 

For each terminal a, and given energy e we define outgoing waves w^ n (e, r) and incoming 
waves w~ n (e, r), where n is the mode number (see Fig.l). In case of a strong magnetic field 

are ES's and n is the LL number. The waves w ± are normalized to carry a unit flux over 
the cross section of the terminal. Choosing the gauge A x = —By, A y = for a given terminal, 
where x and y are the longitudinal and transverse coordinates in this terminal, one can repre- 
sent the waves as follows: w£(e, r) = exp[i/c+(e)x]0+(e, y), w~(e,r) = exp[ik~(e)x](j)~(e, y). 

In what follows we use "hat" to indicate the magnetic field inversion. It means for 
example that if w£ is an outgoing wave for the field B (i.e. an outgoing ES for LL n) then 
w+ is an outgoing wave for field —B (i.e. an outgoing ES for the same LL near the opposite 
boundary). It is easy to check that w^ n (e, r)* = w^ n (e, r) or equivalently fc^(e) = — k^(e) 
and (j)t(e,y) = i%(e,y)*. 

Different functions (j){y) corresponding to the same wave vector k are eigenfunctions of 
the same Hamiltonian and are orthogonal. This is not the case when two functions 4>i(y) 
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and 4>2{y) correspond to the same energy e, but to different wave vectors k\ and In this 



case the "orthogonality" relations are [18 



/ dyfaMih - -A x ) + (k 2 - -A x )\ = 0, (Al) 

and 

/ dyfaUih - -A x ) - (h + -A,)] = 0. (A2) 

For a given energy e the incoming field in terminal a is a superposition of incoming waves 
X) n a m (e)w~ n (e,r), while the outgoing field in terminal (3 is a superposition of outgoing 
waves J2m ^/3m(^)wt m (e, r )- The scattering matrix connects the amplitudes of the incoming 
and outgoing waves 

m,atn (e)a an (e). (A3) 

an 

The scattering matrix is unitary due to flux conservation 

(A4) 

/3m 

and due to time reversal 

(e) = S an ^ m (e). (A5) 

A scattering state x an (e,r) is defined as a solution of the Schroedinger equation with 
energy e excited by an incoming wave w~ n (e, r). Complex conjugate scattering states are so- 
lutions of the Schroedinger equation with inverted magnetic field. Comparing the behaviour 
of x* an d X a t infinity one finds that 

Xan(e, r)* = J2 S f3rn,an(e)* X/3m(e, r) (A6) 
/3m 



and also 



Xan(e, r) =J2 S/3m,an(e) X/3m(t, r)* . (A7) 
ffm 



13 



The Green function is defined by the equation 

[H(r)-e]G e (r,r') = -5(r-r'), (A8) 

with the Hamiltonian given by Eq. ([]]). The boundary conditions are G e (r,r') = when r 
is at the boundary of the NS, and these correspond to outgoing waves when r approaches 
infinity in some of its terminals. The Green theorem in case B ^ is as follows [Tj| 



/ (r)(uHv — vH*u) = (j> dl n 



6 6 

u(—iV A)v - v(-iV + -A)u 

c c 



(A9) 



where n is the unit normal vector directed outside the NS, dl is an element of the boundary. 
Using this theorem and Eq. (|A2|) one can prove the symmetry G t (r', r) = G e (r, r'). 
When r approaches infinity in terminal f3 

Xan (e, r) | r ^oo/3 = S a pw~ n (e, r ) + 53 Sp (e,r), (A10) 

m 

and 

G^r.r)!™^ = -iY, w pme( r )Xpme(r'). (All) 

m 

The first equation follows from the definition of the scattering states and the scattering 
matrix, while the second equation can be obtained as follows. From the explicit expression 
of the Green function for a waveguide in a magnetic field given in ]20| one can see that a unit 
source — S(r — r') at r' — > oo/3 excites an incoming field —iJ2m w pm( e ^ r ')* w pm( e ^ r )- Since 
each wave Wg m (e, r) excites a state Xpm{ e i r ) we fi n d that 

G e (r,r')|r'^oo/3 = -i^Wp m (e,r')*xpm(e,r)- (A12) 

m 

Using the symmetry of G, and the relation between w~ and w + , we find the relation given 
above. 

A useful function is defined as follows 

G e (ry)-G e (r,r')* = -ig e (r,r'). (A13) 

This function can be presented in terms of the scattering states 
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(A14) 



an 



Obviously g e (r, r) is the local density of states. Inverting B one finds g e (r,r') = g e (r,r')*. 
For B = the function g e (r, r') is real. 

Let the confining potential U(r) be subjected to some variation 5U(r). The variation of 
the scattering states Sxam(^, r) contains only outgoing waves and can be found from the first 
Born approximation using the retarded Green function corresponding to the potential U(r). 
The asymptotic behaviour of 5xam(t, r) at r — ► oo/3 can then be found using Eq. (|All|) . As 
a result the variation of the scattering matrix is 
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(A15) 



SU(r) 
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FIGURES 
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FIG. 1. The quantum dot and the nanostructure. w^ n are waves emitted and absorbed by 
terminal a (see text). 
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FIG. 2. A two-terminal device. Scattering states emitted from the terminals are labeled by 
the terminals' corresponding numbers. The interaction area is shown. 
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